SOLVABLE MODELS FOR THE SCHRODINGER 
OPERATORS WITH J'-LIKE POTENTIALS 



YURIY D. GOLOVATY, STEPAN S. MAN'KO 

Abstract. We turn back to the well known problem of interpretation 
of the Schrodinger operator with the pseudopotential a5'(x), where 
5{x) is the Dirac function and a G R. We show that the problem in its 
conventional formulation contains hidden parameters and the choice of 
the proper selfadjoint operator is ambiguously determined. We study 
the asymptotic behavior of spectra and eigenvectors of the Hamilto- 
nians with increasing smooth potentials perturbed by short-range po- 
tentials a.e~^'^{e~^x). Appropriate solvable models are constructed 
and the corresponding approximation theorems are proved. We intro- 
duce the concepts of the resonance set and the coupling function 
6^1 : — > R, which are spectral characteristics of the potential 'P. The 
selfadjoint operators in the solvable models are determined by means 
of the resonance set and the coupling function. 



Introduction 

In this paper we investigate the one-dimensional Schrodinger operators 
with the coefficients having singular support on a disjoint set of points. 
Such operators are named operators with singular point interactions and 
are closely related to exactly solvable problems in quantum mechanics, 
atomic physics and acoustics. Such point interaction models are solvable 
in the sense that the resolvents of corresponding operators can be given 
explicitly. As a consequence the spectrum and scattering quantities can 
also be determined explicitly. 

The main question here is how to interpret the differential operators with 
coefficients equal to the generalized functions since space of distributions 
I?'(M") is not an algebra with respect to the "pointwise" multiplication. 
Generally the problem of choosing a single Hamiltonian corresponding to 
a particular singular point interaction is complicated enough. The theory 
of one-dimensional Schrodinger operators with singular potentials has been 
summarized in research monographs [jj i3j . 

According to the postulates of quantum mechanics every physically mea- 
surable quantity is specified by a selfadjoint operator in a Hilbert space. 
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Starting with the pioneering work of F. Berezin and L. Faddeev [8], the sin- 
gular point interactions have been studied with the help of the selfadjoint 
extensions theory for symmetric operators. As a rule, the symmetric oper- 
ators in physical models are generated by differential operators. This leads 
to the question how to describe all selfadjoint extensions in terms of the 
boundary value conditions. In the case of ordinary differential operators 
the problem has been solved by M. Krein [7} I15j. 

For partial differential operators the problem becomes more complicated 
by the infinite deficiency indices of symmetric operators. In that case new 
approach to the theory of selfadjoint extensions is based on the concept of 
boundary triples [23\ I24j and the corresponding Weyl functions. In [161 117j 
the concept of Weyl function was generalized to an arbitrary symmetric 
operator with infinite deficiency indices. It is well known that the Weyl 
function can be used to investigate the spectrum of a selfadjoint extension 
and the scattering matrix [9|. It is worth to note that the concept of 
boundary triples is useful in the investigation of the point interaction mo- 
dels [SSIIMIETIEH]- 

The singular finite rank perturbations have been studied in papers [21 
[39| HOl H6] and monograph [IT] . The Schrodinger operators with different 
singular potentials have been investigated in [311 dSl [HI [531 15^155] . 

In some instances the set of all selfadjoint extensions of a symmetric op- 
erators is rich enough. Therefore the harder question comes: how to choose 
an extension that is best suited to our physical model. For some models 
the proper operator can not be chosen within the selfadjoint extensions 
theory, because the models contain hidden parameters. If we replace the 
singular potential with a sequence of short-range smooth potentials, then 
the operator obtained in the zero-range limit, as often happens, can depend 
on the type of regularization, that is to say, the operator is governed by the 
profile of squeezed potentials. This profile is a hidden parameter and plays 
a crucial role in the choice of a selfadjoint extension corresponding to the 
physical model under consideration. 

Differential operators with singular coefficients can be studied within the 
framework of the contemporary theories such as the Colombeau algebras 
[12j or the Egorov generalized function algebra [28]. These theories allow 
us not only to multiply the generalized functions, but they also contain rich 
sets of "(5-functions with fixed profiles" as different elements of the algebra. 
This approach to study of solvable models has been used in [5l [6] . 



1. Problem Statement and Main Results 

1.1. Problem statement. In this work we turn back to the well-known 
problem how to define the one-dimensional Schrodinger operator with the 
(^'-potential. New attempt will be made to answer this question. 
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Let us consider the formal Hamiltonian 

Ha = — z-^ + U(x) + a6'(x), xGR, 

where C/ is a real smooth function, 6' is the first derivative of the Dirac 
delta function, and a is a real constant. If product 5'{x)y{x) is treated as 
y{Q)5'{x) — y'{0)6{x), then for a 7^ there exist no solutions in P'(R) to 
equation Hay = Ay, except for a trivial one. In order for operator Ha to 
be assigned a meaning, firstly we must describe all selfadjoint extensions 
in L2(M) of the symmetric operator 

L = + U{x), V{L) = {/ G Co-(M) : /(O) = /'(O) = 0}. 

Next, one of the extensions must be prescribed for the Schrodinger op- 
erators with potential U{x) + a5'{x). At present, in both physics- and 
mathematics-oriented literature there is no a consensus regarding the pro- 
per choice of such operator, since L has a rich set of selfadjoint extensions 
(see Sec. m for the historical remarks). 

Let us consider more realistic model, namely the family of Hamiltoni- 
ans with smooth short-range potentials, which approximate the singular 
potential a5'{x). Let 7^e(a,^') be the closure in L2(M) of the essentially 
selfadjoint operator [ini p. 50] 

F,(a, m) = -^ + U{x) + ^*(e-^x), V{He{a, ^)) = Co~(M). 

Here e is a small positive parameter. We call ^ G C^(M) the profile of a 
local perturbation, and a the coupling constant. It is clear that for some 
profiles sequence e~^^(e~"^x) converges to 5'{x) in !D'(M). We will denote 
by V the set of real functions ^' G Co°(M) such that supp^' = [-1, 1]. Let 
£{L) denote the set of all selfadjoint extensions of L. Suppose also that 
potential U increases as |x| 00. It follows that the spectrum of TCe{ct, ^) 
is discrete. 

Our main goal is to construct the map R x V — > ^{1^) that assigns 
a selfadjoint extension TC{a,'^) of operator L to each pair (a,^). The 
choice of an operator is determined by the proximity of the energy levels 
and the pure states for the Hamiltonians with smooth and singular po- 
tentials respectively. We explore the asymptotic behavior of eigenvalues 
and eigenfunctions as e — > 0. The asymptotics provides the limit operator 
W(a,^'). 

1.2. Structure of Paper. The paper is organized as follows. In SecEJ we 
derive qualitative properties of the discrete spectrum of perturbed operators 
We (a, ^). We show that all eigenvalues are continuous functions on e that 
are bounded from above. Generally speaking, the spectrum of this family is 
not bounded from below. For some profiles ^ and a ^ there exists a finite 
number of eigenvalues which go to —00 as e ^ 0. By the bounded spectrum 



4 



YURIY D. GOLOVATY, STEPAN S. MAN'KO 



of Tieioj^if) we mean the set of such eigenvalues that remain bounded as 
e ^ 0. 

The leading terms of asymptotic expansions for the eigenvalues of the 
bounded spectrum and the limit operators TC{a, ^) are formally constructed 
in Sec. [3l We introduce two spectral characteristics of profile namely 
the resonance set Sij, that is the spectrum of the Sturm-Liouville prob- 
lem —w" + a^w = 0, w'{—l) = 0, w'{l) = on interval (—1,1) with 
respect to spectral parameter a, and the coupling function O-q/ : T^-q, M. 
We set 04f(a) = u)q,(1) (tt;Q,(— 1)) , where Wa is an eigenfunction corre- 
sponding to eigenvalue a £ S^r. In the case, when the coupling constant 
doesn't belong to the resonant set, 7i{a, ^) is just the direct sum of the 
Schrodinger operators with potential U on half-axes subject to the Dirichlet 
boundary condition at the origin. In the resonant case, when a £ Sij,, op- 
erator H{a, is the closure in L2(M) of the essentially selfadjoint operator 
Hia,"^) = + U with domain P(i7(a,^)) = {/ S /(+0) = 

e^{a)f{-0), 4(«)/'(+o) = /'(-o)}. 

In Sec. m we provide a detail discussion of the previous results dealing 
with the J'-interactions and the ^'-potentials. We give two examples of 
exactly solvable models with piecewise constant J'-like potentials for which 
the resonance set and the coupling function are explicitly computed. We 
solve the scattering problem for this J'-like potential. The limit value of the 
transmission coefficient can be also explicitly computed in terms of coupling 
function 0^. 

The remainders of asymptotics for the eigenvalues and eigenfunctions of 
7ie{a., ^) are constructed in Sec. [5l because we are in need of more precise 
asymptotics, in order to proof the approximation theorems. 

Section [6] contains the main results of this paper. We prove in this sec- 
tion that the eigenvalues of ^^(q;, ^) that belong to the bounded spectrum 
converge to eigenvalues of T-L{a, ^). The convergence in L2(M) of the cor- 
responding eigenfunctions is also established. 

2. Spectrum of He{a,^) and Auxiliary Results 

Firstly we describe the set £{L) of all selfadjoint extensions of minimal 
operator L. The adjoint operator L* = —-^ + U{x) is defined on the 
domain 

V{L*) ={v£ \ 0) : -v" + Uv£ LaW}. 

Lemma 2.1. Every element of £{L) coincides with operator L* , restricted 
to the set of functions, satisfying the boundary conditions at the origin of 
one of the types 

(2.1) h^v'{-Q) = /i^u(-O), h+v'{+Q) = /i^f(+0) 
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with the parameters {hf,h^) from the projective space P^; 

(-) t'^^) - <:'(-o)) • ^-"(2 z 

where G [— §, Cki G M and C11C22 — C12C21 = 1. 

This lemma has been proved in |1H [50] for {7 = 0. Obviously, the 
smooth potential does not affect on the boundary conditions at the origin. 
Selfadjoint operators described by boundary conditions (|2.ip will be named 
separated extensions. These operators are equal to the orthogonal sum 
of two selfadjoint operators defined on the half-axes. Selfadjoint opera- 
tors of the second type will be named connected extensions, because these 
conditions connect the boundary values of the function on the left and right 
half-axes. 

Let us write ^'e(x) = e'^^{e^^x) and set mk{f) = /^^^/(O ^C- 

Lemma 2.2. Assume ^ £ V and c is a nonzero constant. Sequence 
converges to c6'{x) as e ^ in the sense of distributions if and only if 

(2.3) mo(^') = 0, mi(^')/0. 

In addition, c = —mi{^). 

Proof. Given ip G C^(IR), we have 

£ 1 

{^s,^)=e-^ [ ^{e-^xMx)dx = e-^ f ^{e)p{eOdi 



1 

^!{i){^{Q)+e^'m + 0{e^))d^ 

-1 



= £- Vo(^')93(0) + mi{^)^'{^) + 0{e) 

as e — > 0. Hence sequence {^s-,'^) converges for all if G Co°(M) to a finite 
limit if and only if mo(^) = 0. Moreover if mi{^) 7^ 0, then the limit is 
nontrivial and the limit functional is defined by ^0 = —'mi{^)5'{x), since 
(^-0,^) =mi(1')(/p'(0). □ 

The first condition in (j2.3p implies that function ^ must change a sign on 
[— 1, 1]. The second one involves an asymmetric property of this function: 
if ^{Cj satisfies (12. Sp . then — a) cannot be an even function for any a. 
We introduce the set of 5' -like potentials 

Vq = {^ £ V: mo(^') = 0, mx{^) = -1}. 

Recall that U{x) — > +00 as |x| +00. Then for all e > the spectrum 
of 'He(a,^') is real, discrete and simple. Let {A|(a,^')}^]^ be the eigenva- 
lues of 'He(a,^') enumerated in increasing. Suppose that a, ^')}^;^ 
is the orthonormal in -L2(M) system of eigenfunctions. 
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Theorem 2.1. For each pairs (a,^) G M x P eigenvalues A|(a, are 
continuous functions with respect to £ £ (0,1). Moreover, the eigenvalues 
are bounded from above as e ^ 0. In the case of the 6' -like potential, the 
spectrum of 7iir{a,^) is unbounded from below as e — > 0; if a ^ and 
^ G Vq, then AKaj^*) < —76"^ for some positive constant 7. There is 
only finite number N{a, ^) of eigenvalues that converge to —00 as e ^ 0. 

Proof. Fix a and ^. The quadratic form 

as[u]= [ {\u\'^ + {U + a^e)\u\'^)dx, u£C^{R) 
Jr 

is a continuous function with respect to e G (0, 1). The minimax principle 
A|.(a, = inf sup 

veEk, \\v\\=l 

yields the continuity of eigenvalues with respect to e. Here -E^ is a k- 
dimensional linear subspace of C^(M), and || • || denotes the norm in L2(M). 

We choose subspace E"^ , whose elements vanish in a neighborhood of the 
origin. Then, we have 

(2.4) A|(a,^')< sup ae[v]. 

v&El, \\v\\=l 

For all sufficiently small e the restriction of to the finite-dimensional 
space -E^ does not depend on e. This yields the boundedness of the eigen- 
values from above. 

From now on, ^' is a (5'-like profile. It is well known [51, p. 338] that the 
Schrodinger operator — ^-|-a^'(^) on the line has a negative eigenvalue for 
all a 7^ if and only if ?7io(^') = 0. Hence there exists function u, ||ii|| = 1, 
such that 

^ = j + a-^{^)u^) d( < 0. 

R 
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Let C £ Ccr'(l^) be the positive cutoff function such that suppC = [—2,2] 
and (^(x) = 1 for x € [—1,1]. For abbreviation, we write Q and in- 
stead of C{e^) and C'i^C) respectively. Let us consider sequence Ue{x) = 
e-^/2^(x)tt(e-^x). Then 



u, 



|2 



e I 



i J C\x)u\l)dx = J C^iC)u\OdC^ J u\OdC = l 



R R a 

as e — > 0. Next, we conclude from the minimax principle that 

\ / X ■ r r 1 ^ ae[Ue] 

^GC7o°°, ^ - \\Ue\\ 
\\v\\=l 

hence that 

e^lluell Ai(e,a) < e'^ae[ue] 

e I {{C{x)u{e-^x)y^ + {U + a^e){Cix)uie-^x)f)dx 



+ I C'eU^di + e I U{x)C^{x)u^{e-^x)dx. 



K -2 

It is easily seen that the first integral converges to negative number /i as 
e — > 0, and other terms go to zero. Thus for e sufficiently small the estimate 
e^lliiell Ai(e,Q;) < ///2 holds, and so Ai(e, a) < —76"^ for some 7 > 0. 

Let N~{a, ^) be the number of negative eigenvalues of 7ie{oi, Obvi- 
ously N[a^^) < N^{a^'^) for e > sufficiently small. In the case of the 
continuous potential the inequality 

(2.5) N-{a,-^)<l + j \x\\U'{x)\dx + \a\ j |x| 

R R 

holds [lOi p. 97], where f~{x) = min{/(x),0} is the negative part of /. 
Since U increases as x — > 00, U~ is a function with compact support, and 
consequently the first integral in (12. 5|) is defined. The support of lies 
in [— e,e], we thus deduce 

£ 1 

/N|*,W|...s-/HI*-(f)M.^/|«||*-K)|.e. 
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In view of ()2.5p and the formula above N(a, ^) is uniformly bounded with 
respect to e: 

N{a, < N-{a, < ci{U) + C2(^')|a| 
with constants ci{U), C2(^') being positive. □ 

Therefore the spectrum of ^{^{a, consists of two parts: {A|(q;, 
is the set of eigenvalues that tend to -co as e ^ 0, and {A|(a, ^')}^^_,_]^ 
is the set of bounded eigenvalues, which we call the bounded spectrum 
of 'He(a,^'). Notice that eigenfunctions yf (x; a, ^'), . . . , y|^r(3;; a, have 
strongly oscillatory character as e — > in an neighborhood of the origin 
and vanish exponentially outside. They converge to zero in L2(K) weakly. 



3. asymptotics of bounded spectrum of we(a,^'): leading 

Terms 

Let us consider the eigenvalue problem 

(3.1) -y'l+{U{x)+ae-H{e-^x))ye = \'ye, Ve^^L^m. 

Fix eigenvalue A|(a, ^) with number k > N{a, We denote it briefly by 
A^. Let ys be the corresponding eigenfunction. We postulate asymptotic 
expansions for the eigenvalue 

(3.2) X' X + eXi+e^X2 + --- , 

and two-scale expansions for the eigenfunction 

(3.3) ys{x) v{x) + evi{x) + e'^V2{x) + ■ ■ ■ for |a;| > e, 

(3.4) ydx) ~ w{e~'^x) + ewi{e~'^x) + e'^W2{e~^x) + • • • for |x| < e. 

Here v, Vi are defined for x G M \ {0} and belong to L2(M). Set ^ = e~^x. 
Functions w, Wi are defined for ^ e [—1,1]. We also assume that v is 
deferent from zero. Series (13. 3p . (j3.4p satisfy the coupling conditions 



(3-5) [y.L=±. = o, [y;],=±, = o 

at points x = ite, where [/]a;=a = /(a + 0) — /(a — 0) is the jump of / at 
point a. This type of combined asymptotic expansions has been used in 
particular in [T8]-p0]. 

Upon substituting (j3.2p - (|3.4p into equation (j3.ip . we obtain 



(3.6) -v" + U{x)v = Xv, xGM\{0}, 

(3.7) -v'i + U{x)vi = Xvi + Xiv, xEM\{0}, 
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and also 

(3.8) -w" + a^{Ow = 0, 

(3.9) -w'l + a^{Owi = 0, ee(-l,l), 

(3.10) -W2+a^{0w2 = Xw-U{0)w, ^g(-1,1). 

Next, substituting series ()3.3p . (|3.4p into the coupling conditions yields 

v{±e) +evi{±e) H ~ w{±l) +ewi{±l) H , 

v'{±e) +ev[{±e) + e-^w'{±l) + w[{±l) + £^2(^1) + • • • ■ 

We can now expand u(ibe), v'{zize) into the formal Taylor series about 
X = ±0. Then 

(3.11) v{-0) = w{-l), v{+0)=w{+l), 

(3.12) w'{-l) = 0, w'{l) = 0, 

(3.13) v'{-0) = w[{-l), v'{+0) = w[il), 

(3.14) vi{-0)-v'{-0)=wi{-l), vi{+0) + v'{+0) = wi{l), 

(3.15) v[{-0) - v'\-0) = w'^i-l), v[{+0) + v"{+0) = w'^il). 

It follows that V satisfies equation ()3.6p on each half-axis and is a 
solution to the problem 

(3.16) -w" + a^w = 0, ^e(-l,l), w'{-l) = 0, w'{l) =0. 

Moreover both these functions satisfy coupling conditions (|3.1ip . Problem 
(|3.16p is decisive in our next considerations, because it contains information 
about a kind of the singular perturbation, namely the profile ^ and the 
coupling constant a. The first and primary question is whether there exists 
its nontrivial solution. 

3.1. Resonant Set and Coupling Function. Problem (j3.16p can be 

viewed as a spectral problem with spectral parameter a. Since function ^ 

is sign-changing, it is natural to introduce a space with an indefinite metric. 

1 

Let /C be the weight L2-space with the scalar product (/,<?) = / |^'| fgd£,. 

-1 

1 

Let us introduce in /C the indefinite metric [f,g] = J ^f'gdS,. Then the 

-1 

pair (/C, [•,•]) is called a Krein space (see [1] for the original definition). In 
the Krein space there exists the fundamental symmetry Jf = sgn ^ ■ f such 
that [f,g] = {Jf,g) for ah /,£/ G /C. 

Let T be a closed densely defined operator in /C. The J-adjoint operator 
of T is defined by the relation [Tx,y] = [x,T^*^y] for x G 2?(T) on the 
set of all y G /C such that the mapping x ^ [Tx, y] is a continuous linear 
functional on T>(T). The operator T is called J-selfadjoint if T = T^*!. The 
operator T is called J-nonnegative if [Ta;,^] > for x G V[T). 
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Lemma 3.1. Let T be a J -self adjoint and J-nonnegative operator with 
a nonempty resolvent set. Spectrum cr{T) is real, moreover, the residual 
spectrum is empty. If X ^ is an eigenvalue ofT, then it is simple (e.i., 
ker(r — A) = ker(r — A)^j. If Q is an eigenvalue ofT, then its Riesz index 
is equal or less than 2, e.i., kerT^ = kerT^. Generally, is not a simple 
eigenvalue [4, p. 138]. 

We introduce in /C the operator 7^ = —-^j^^ with the domain 

P(T*) = {/ e /c I / G wii-l, 1), ^^V" e /C, /'(-I) = 0, /'(I) = 0}. 

Spectral equation T^w = —aw is associated with (13.16j) . Differential equa- 
tions with the indefinite weights have been widely studied [141 [32| \33 \ I49 j . 
For the sake of completeness we describe the main properties of Tvj-. 

Theorem 3.1. For each ^ £ V operator Ti^ is J -self adjoint and J-nonne- 
gative. 



Proof. Given / E D(Tq,), we have 



1 



[T^f,g] = - 1 f"gd^ = f{l)g'{l)-f{-l)g'{-l)- J f g" d^. 



-1 



Then identity [T^f,g] = [/.T^^] holds if 5 G V{T^) and T^^ = -m~^g". 
Hence 7^ is J-selfadjoint. Next, for all / G V{T^^) 

1 1 

[r^fJ] = - j f'JdC = J |/f dC > 0, 

-1 -1 

so that 7^ is a J-nonnegative operator. □ 

Theorem 3.2. (i) Given ^ £ V, the spectrum of Tq, is real and 

discrete. All nonzero eigenvalues are simple, //mo(^) / 0, then 
is a simple eigenvalue. In the general case, kerT^ 7^ kerT^. 
(ii) If^ G Vq, then the spectrum ofTq; has two accumulation points 
—00 and +00. Moreover, kerT^ 7^ kerT^. 

Proof. To start with, show that the resolvent set of 7^ is nonempty. Each 
solution of the homogenous problem 

(3.17) / + i^r^ = 0, ^G(-l,l), g'{-l)=0, g\l)=0, 

satisfies equality 

1 1 
\g'fd^,-i I ^'|g|2d^ = o. 



-1 -1 

Since function ^ is real-valued, we deduce that g is a constant function. 
Obviously, the zero function can only be a solution of ()3.17p . Then the 
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nonhomogeneous problem g" + i^g = f, g'{—l) = 0, ^'(l) = admits a 
unique solution for each f £ JC. Hence i G p{T^) and g = R{i,Tiii)f, where 
R{X,Tq/) is the resolvent of T^. 

Therefore spectrum u{T^) is real by Lemma l3. II and Theorem 13.11 The 
compactness of R{i,Tq,) follows from the sequence of embeddings 



As a consequence we have o"(7^) = (Tp(7^). 

Each eigenvalue of 7^ is simple. Indeed, if there exist two linear indepen- 
dent eigenfunctions <p and tp corresponding to eigenvalue a, then function 
1^(0 = ~ '/^(l)V'(O is nonzero since both values 97(1), ■0(1) are 

different from 0. But this is impossible, because = as a unique solution 
to the Cauchy problem w" + a'ifw = 0, (1) = 0, = 0. 

Assume now £ Vq- Since potential ^' changes a sign, the spectrum 
(t(7^) is unbounded in both directions [H] (see Remark 13. ip . Moreover 
invariant subspace ker is generated by eigenvector w = 1 and root vector 



which is a solution to the problem = — ^'(^), w^— 1) = 0, w'^{l) = 0. 
This problem admits a solution if and only if ?7io(^') = 0. In view of 
Lemma [3. II there are no linear independent root vectors other than w^:. □ 

Remark 3.1. The spectrum of Tg, is a symmetric set with respect to the 
origin for each odd profile ^: if a is an eigenvalue with eigenfunction w(^), 
then —a is also an eigenvalue with eigenfunction w{—^). 

We introduce the set = {a G M: — a G a{T^)}. We call E^r the reso- 
nant set of profile ^. Problem ()3.16p has a nontrivial solutions if and only 
if the coupling constant a belongs to Suppose Wa is an eigenfunction 
of (|3.16p with eigenvalue a S S^,. Thus values Wq,(— 1) and Wq,(1) are 
nonzero and so the quotient 



is correctly defined. Quantity 6\i/{a) does not depend on the choice of an 
eigenfunction and can be considered as function 6qj : M defined at 

points of the resonant set. We call 9\i, the coupling function of ^. Hence 
each profile ^ generates the resonant set and the coupling function. 



V{T^) C W^i-l, 1) ^ L2(-l, 1) C /C. 




1 



3.2. Limit Operator TC{a, We now continue to construct the asymp- 
totics. We distinguish two different cases. First assume a does not belong 
to resonant set . Then problem (j3.16p has a trivial solution w = only. 
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and so from coupling conditions ()3.1ip we deduce v{—0) = ^(+0) = 0. 
Recalling (|3.6p . we obtain 



(3.18) -v" + Uv = Xv, x£R\{0}, u(0)=0, ^£^2(1^). 

Function v is nonzero by our assumption, thus A is an eigenvalue of the 
problem ()3.18p . Therefore we can introduce limit operator 7i{a,^), which 
is the direct sum S- © S+ of the Schrodinger operators acting on half-axes. 
Operators S- and S+ are associated with problems 



(3.19) 



-v" + Uv = Xv, xGM_, j -v" + Uv = Xv, x G M+, 

v{0) = 0, ve L2{R-), \ v{0) = 0, ve L2(M+) 



respectively, where M_ and M+ are the negative and positive half-axes. 
Spectra of 5_ and 5"+ are real, discrete and simple. Obviously, 

a{n{a,^)) = a{S-)Ua{S+). 

Let us now suppose that a belongs to resonant set Ti^j. We set w = aWa, 
where Wa is the corresponding eigenfunction of ()3.16p . and a / 0. We 
conclude from (|3lT]) that v{-0) = aWa{-l), v{+0) = aWa(l), hence that 

(3.20) v{+0) = eii,{a)v{-0). 

According to (13. 9p . (j3.13p the next term wi of series (13. 4j) can be found by 
solving the problem 

(3.21) -w'l + a^wi = 0,^e (-1,1), w[{-l) = v'{-0),w[{l) =v'{+0). 

Because a is an eigenvalue of (|3.16p . the problem admits a solution if and 
only if Wail)v' (+0) = Wa(— l)f'(— 0). This solvability condition can be 
derived multiplying the equation by Wa and integrating twice by parts. It 
may be written in the form 

(3.22) eq,{a)v'{+0) = v'{-0). 

Combining (j3.6p . (|3.20p and (|3.22p we deduce that v must be an eigenfunc- 
tion of the problem 



(3.23) 



-v" + Uv = Xv, x£R\{0}, 

v{+0) - 9ii,{a)v{-0) = 0, 9^{a)v'{+0) - v'{-0) 



Thus the connected selfadjoint extension of operator L 

n{a,^) = -—r + Uix), 
(3-24) V{'H{a,^)) = {feV{L*):f{+0) 

= e^ia)f{-0), 9^{a)f' i+0) = /'(-O)} 

is associated with (|3.23p . 
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Remark 3.2. In the case a G Sij, we can also assume that a = 0, that 
is, w = 0. Then t'(O) = 0, and so v is an eigenfunction of S'_ © 5"+. But 
condition (j3.22p has nevertheless to hold true; for otherwise we could not 
solve (I3.2ip . If V is an eigenfunction corresponding to a simple eigenvalue 
of S- © 5"+, then at least one of values 0) and i;'(+0) equals 0, since v 
vanishes on at least one of half-axes. Then condition (j3.22p yields f '(— 0) = 
f'(-|-0) = 0, which is impossible. For a double eigenvalue of S- © 5+, there 
exists a unique line in the 2-dimensional eigenspace, for which condition 
(j3.22p can be satisfied. Vectors of this line are also eigenfunctions of (|3.23p . 
since condition (j3.20p is trivially valid. Therefore in the resonant case the 
leading term w of (j3.4p is zero if and only if A € a(Ti.{a, ^)) n a{S- © 
This intersection consists of all double eigenvalues of S- (B S^. Note that 
for operator 7i{a,^) these eigenvalues are simple. 

Remark 3.3. For a = operators TCeiO, ^) don't depend on e and coin- 
cide with the unperturbed Schrodinger operator with potential U. On the 
other hand, the resonant set always contains 0, and so operator 0.(0,^) 
must be defined as in (j3.24p . Obviously 6q,{0) = 1, because the constant 
eigenfunction corresponds to the zero eigenvalue of 7^. It follows that 
TisiO, = n{0, ) for ah e > 0. 

After all, for each profile ^ we have defined the family of self adjoint 
operators {7i{a, ^)}aeR- If coupling constant a doesn't belong to resonant 
set Sij,, then operator Ti.{a, is a direct sum of the Schrodinger operators 
acting on half-axes. The corresponding eigenfunctions describe states of a 
quantum system, when a particle can be found with probability 1 on one 
of the half-axes. In this case we obtain what will be referred to as the 
closed 5' -barrier. If a G T,xj,, then operator Ti.{a,^) given by (|3.24p is the 
connected selfadjoint extension with matrix 



in (12. 2p . For resonant coupling constants a particle can permeate through 
the barrier with some probability (the case of an open 5' -barrier). 



4.1. Historical Remarks. The formal Hamiltonians with the Dirac func- 
tions and their derivatives in potentials were studied since the eighties of 
last century. We refer, for example, to research monograph [1] (the first 
edition, 1988) and papers [50l [52l [26] . Incidentally, in solvable models one 
distinguishes two different phenomena dealing with the (^'-function: the 
S' -interaction and the point dipole interaction (the 5' -potential) . 

The first result about the (5'-interaction is due to S. Albeverio, F. Geszte- 
sy, R. H0egh-Krohn and H. Holden ^j. The standard definition is the 



(3.25) 




4. What Is the (5'-Potential? 
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one-parameter family of selfadjoint extensions 

A =- — 

(^•1) V{Ap) = {/ G VF|(R \ {0}) : /'(-O) = /'(+0), 

/(+0)-/(-0)=/?/'(0)}. 

P. Seba |52j showed that these selfadjoint extensions correspond to the 
heuristic operator 

-^^+fiW{x)){6'{x)\, 
where the rank one perturbation |/) (/| is defined by 

{\f){fW)){x)=f{x) I f{yMy)dy. 

R 

In the same paper P. Seba also studied formal Hamiltonian —-^ + a6'{x), 
where 5'{x) was considered as the limit in the topology of P'(M) of linear 
combination 

^{6{x + e) - 6{x - e)) . 

The limit operator is referred to as the (5-interaction dipole Hamiltonian. 
The author proved that this Hamiltonian is the direct sum S- © 5"+ with 
U = (see (j3.19p ). For instance, in Theorem 4 one proves that the family 
of the Schrodinger operators with smooth potentials 

(4.2) s.(V) = -^^j,v(^). ..R 

converges to S- © 5"+ in the norm resolvent sense, provided potential V G 
C^(M) fulfils condition mo(T^) = 0. From viewpoint of the scattering 
theory this means that the 6' is a totally reflecting wall. 

But these conclusions are at variance with the results presented in [13^ 
[561 [581 [59] , where the scattering properties of 6' -like potentials constructed 
in the form of squeezed rectangles were studied in the zero-range limit 
(a simple such problem is treated in Subsection 14. Sp . The resonances in 
the transmission probability were firstly described in [13]. The authors 
obtained a discrete set of resonant coupling constants On for which there 
is a partial transmission decreasing as a„ is larger. These values are roots 
of a transcendent equation, which depends on the regularization, that is to 
say, the profile of piecewise constant potentials. 
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It is worth to point out preprint |42j . P. Kurasov and N. Elander pro- 
posed to define the Hamiltonian with the (^'-potential as the second deriv- 
ative operator in L2(M.) subject to the couphng conditions 

/(+0)-/(-0) = |(/(+0) + /(-0)), 

(4.3) 

/'(+0) - f'i-0) = -|(/'(+0) + /'(-O)). 

This definition is based on the extension of Dirac's functions to the case of 
discontinuous test functions [43] : 

{6^-Hx),^{x)) = ^(c^W(+0) + ^W(-0)). 

The same definition was introduced by L. Nizhnik [l8], who studied the 
Schrodinger operator with the (5'-potential in Sobolev space (K \ 0). 

Alternative approaches to the problem can be found in [251 1301 135j . In 
[31 p. 339] we read that the relation between the set of point interactions 
(selfadjoint extensions from Lemma [2. ip and the set of interactions defined 
by the (5'-potential can not be established without additional assumptions 
on the symmetry properties of the interaction. Different assumptions can 
lead to different one dimensional families of point interactions correspond- 
ing to the formal expression —-^ + a6'{x). It favours the view that this 
problem contains hidden parameters. 

4.2. Solvable Model for Hamiltonian 

with Potential ae^^'^le^^x). The formal asymptotic results obtained in 
the previous section point to the following fact: the solvable model that is 
appropriate to operator —-^ + a6'{x) is not specific for the given point 
interaction, it depends on the way in which the derivative of Dirac's func- 
tion is approximated in the weak topology. In fact, profile ^ is a hidden 
parameter in the conventional formulation of problem how to define the 
(^'-potential. 

However, the solvable model that most closely corresponds to the motion 
of a quantum particle in potential p-^'(e~^x) is uniquely defined. The model 

is described by the family of operators A(a,^) = —-^ with domains 
V{A[a, ^)) = {fe W^{R \ 0) : /(-O) = /(+0) = 0} for a S.,, 
V{A{a, ^)) = {/ G VF|(M \ 0) : /(+0) = e^{a)f{-0), 

0*(a)/'(+O)=/'(-O)}/orQGSv,, 

where and are the resonance set and the coupling function of ^ . 

This solvable model is in good agreement with the results obtained previ- 
ously in [m [Ml [SSI [59] . We actually describe the phenomenon of resonance 
in the transmission probability for arbitrary profiles ^ and obtain the res- 
onance set as a spectral characteristic of the profile. 
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It is worth to compare our definition with that proposed in [l2l US] (see 
4.3p ). in which the couphng matrix 



'2 









2-a 
2+a 



is also diagonal provided \a\ ^ 2. In addition, for two values of a the 
model is described by separated selfadjoint extensions, namely /'(— 0) = 0, 
/(+0) = for a = -2, and /(-O) = 0, /'(+0) = for a = 2. 

From the physical viewpoint smooth potential U has no effect on the 
coupling conditions at the origin. Moreover, the resonance set and the 
coupling function don't depend upon U. We cannot but assume that U = 
in our computations above, because operators TCe{a, lose discreteness of 
the spectrum. We shall give two examples of exactly solvable models with 
piecewise constant (5'-like potentials in order to motivate our definition of 
the point dipole interaction. 



4.3. Transition of Quantum Particle 

Through 5'-Like Potential. Let us consider the equation 

(4.4) -y" + ae~'^'^{e'^x)y = k^y, x G M, 

with ^ being a 5'-like profile, such that ^'(^) = 1 for ^ G (—1,0) and 
^'(^) = —1 for ^ G (0,1), A; > 0. Clearly the results obtained above are 
still true for the piecewise constant potentials. Assume that a = > 0. 
It is desired to find solution y^{x;>c,k) of (j4.4p . such that 7/e(x;x, fc) = 
^ikx _j_ Ji^{^x , k)e~^^^ for x < —£ and ye{x;>i,k) = T^{>i , k)e'^^^ for x > e. 
Quantities |i?e(K, /i:)p, \Tfr{K,k)\^ are called a reflection coefficient and a 
transmission coefficient respectively. They have a probabilistic meaning in 
view of \Rs{>c,k)\^ + \Tfr{x,k)\^ = 1. The same example is treated in [13J. 
We adduce this example to show the direct relationship between scattering 
data and the coupling function introduced above. 

For the profile under consideration resonant set is symmetric with 
respect to the origin, and its positive part consists of the roots of transcen- 
dent equation h{^/a) = 0, where h{x) = x(tanhx — tanx). The coupling 
function is given by 

/A r:\ n f \ COsh ^ -r ^ n n I \ COS 

(4.5) 6l^(a) = ^ if a > 0, 6*^(0) = — ^ ^ if a < 0. 

cos-v/a cosh V— a 



Let ui{^; X, t), U2{^; x, t) be the fundamental system of solutions of 
-n" + (x2^'(^) - r2)u = 0, ee(-l,l), 



such that its elements are smooth functions of x and r. Such solutions 
can be constructed explicitly with trigonometric and hyperbolic functions. 
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Then we can find solution ys{x; x, k) in the form 

f ^tkx _|_ ^-tkx £qj. ^ ^ 

Ce,iiii(f ,x, efc) + C£^2'U'2(f ,x, eA;) for |x| < e, 



Te e"-^"^ for X > e. 



where we intend to select the coefficients to satisfy the C^-smoothness con- 
ditions at points x = ite. A direct computation gives 

^ ' ^ ~ (2ieA; - cos xcosh x + 0{e'^k'^) 

as eA; — > 0. The asymptotic behavior of the transmission coefficient as 
e ^ depends upon whether a = belongs to the resonant set. Namely 

(4.7) |T,(x,fc)|^= .(i + 0(gfc)) forxVS*, 

n'^[>c) cos^ X cosh x 

(4.8) |Te(x,fc)|'= , ^2 •(l + 0(gfc)) forx^GSviy. 

cos^ X cosh X 

Whence it follows that the transmission coefficient goes to zero outside the 
resonant set, and it tends to a positive limit for a G Yi^i. This limit |T(a)p 
is independent of k and can be represented in terms of the coupling function 

4g|(«) 

This formula is still true for negative a. According to (j4.5p . transmission 
probability |T(a)|^ tends to zero as \a\ — > +oo. 

Therefore operators A{a, ^'), defined in the previous subsection, are the 
connected selfadjoint extensions if and only if a E S^. 

4.4. Sturm-Liouville Operator With (5'-Like Potential. Let (o, h) be 

a bounded interval, which contains the origin. We consider the eigenvalue 
problem 

(4.9) -y" + ae-H{e-^x)y = \y, x € {a,b), y{a) = 0, y{b) = 0, 

where ^ is the piecewise constant 5'-like potential as in the previous sub- 
section. The fundamental system of solutions of (j4.9p can be constructed 
explicitly with trigonometric and hyperbolic functions. 

We perform the calculations for positive a only. Let us introduce the 
notation lv = \/A, x = ^/a, /ii(x) = tanhxtanx — 1 and ^(x) = (1 + 
tanh xtan x)(l— tanh xtan x)~^. The characteristic determinant A(e, x; to) 
of (j4.9p . whose roots are eigenfrequencies lOs, admits the asymptotic repre- 
sentation as — > 0: 

(4.10) A(e,x;w) = /i(x) jtan au; tan 6w + etij (tan 6w — tan aw) } 

+ euj hi (x) (tan buj — g{x) tan aoj) + 0{e'^uj'^) 
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with h = h[x) being the characteristic determinant of problem (I3.16P de- 
fined in 14. 3[ 

Suppose that a Then h{>c) is different from 0, and so bounded as 
e — > eigenfrequencies oj^ of (14. 9[) tend to the roots of equation 



tan auj tan ftw = 0. 

The corresponding eigenfunctions converge in C(a, b) and the set of their 
possible limits consists of two families of functions 

'sin^(x - a), x G (a,0) 
0, X G (0,6) ' 



yk,2ix) 



xG(a,0) 
sin^(x - 6), X e (0,6) 



b 

where A; G N. It is easy to check that the limit eigenfrequencies and eigen- 
functions correspond to the direct sum of the second derivative operators 
on intervals (o, 0) and (0, 6) subject to the Dirichlet conditions. The sum 
is similar to operator A(a, ^) for a Sij,, defined in 14.21 

Now assume a G Sij,, and so /i(x) = 0. Remark that at the same time 
hi{x) = tanh^ x — 1 < 0. In this case the eigenfrequencies uj^ go to the 
roots of equation tan buj = g{>c) tanacj with 

. , 1-1- tan^ X cosh^ x n2 / \ r ^ 

g(x) = TT- = ^ — = Oiia) for a G n M+. 

l-tanh^x cos2j< 

The eigenfunctions converge in L2{a,b) as well as uniformly on each of 
intervals (a, 0) and (0, 6) towards functions of the family 



r(u;fc) sinc<jfc(x — a), x G (a, 0) 
^(a) sinu;jfc(x — 6), xG(0,6)' 



where u>k is a root of tan 6a; = ^^(a)tanaa;, and r{u>) = for sin aw ^ 

0, r{oj) = ^ for sin aw = 0. We check at once that w^ and are 
eigenvalues and eigenfunctions of the problem 

(-y" = Xy, X G (a,0)U(0,6), y{a) = y{b) = 0, 
\ y{+0) = 9^{a)y{-0), e^{a)y'{+0) = y'{-0), 

The operator associated with the problem is similar to operator A(a, ^) in 
the resonant case. 

4.5. Open Problem. In Section [3] we have constructed the asymptotics 
in the case of an arbitrary profile belonging to V. The case of 5' -like 
potentials ^' G "Pq is characterized by a special structure of resonance set 
Sij, described in Theorem 13. 2l /'ii ) as well as a behavior of coupling function 
Oxj/. The exactly solvable models above and the computer simulation of 
more complicated models suggest the following 
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Hypothesis. Suppose that potential ^ fulfils conditions mo(^') = 0, mi(^') 
= — 1. Then the coupling function 9\i, possesses the following property: 

o \0^{a)\ > 1 for a £ n M+ and \9\s;{a)\ — > +oo as a ^ +oo, 
o \9ii/{a)\ < 1 for a £ Sij, H M_ and ^0 as a ^ —oo. 

The property has its origins in the structure of eigenfunctions of J- 
selfadjoint and J-nonnegative operator 7^. Because of this, the demon- 
stration could be obtained using the Krein space theory. Note that our 
hypothesis for the behavior of fails for other classes of potentials. There 
exist even functions ^ such that mo(^') 7^ as well as mo(^) = mi(^) = 0. 
But for an even potential we have |0\i,(a)| = 1 for all a € S^. 

The physical interpretation of a coupling function is fairly simple. Sup- 
pose that w is a normalized in L2(M) eigenfunction of (13.231) and Pv{a, b) = 
la '^^ probability of finding a particle in interval (a, b) provided 

that the system is in a pure state v. Then 

which is to say, 9'^ (a) is the marginal ratio of probabilities of finding a par- 
ticle in (0,r) and (— r, 0) respectively. Let us consider potential ^'^{e~^x) 
with the profile being an odd function as shown in Pic. [1] from the left. 
The hypothesis states that the probability to find an electron nearby the 
origin is higher from the side of the deep well than that of the high wall. In 
other words, for positive a a negative drop in the electron density occurs, 
while passing across the junction from the barrier to the well side. If a is 
negative, then positions of the barrier and the well are changed over. 

5. ASYMPTOTICS OF BOUNDED SPECTRUM OF We(a,^'): CORRECTORS 

We need to find next terms of series (|3.2p - (j3.4p in order to justify the 
proximity of energy levels of Hamiltonians 7ie{ot, ^) and 7i{a, 

5.1. Case of Closed ^'-Barrier. If coupling constant a doesn't belong 
to S^, then w = 0. Let A be a simple eigenvalue of 7^(0, = S- ® 5*+ 
with L2(M)-normalized eigenfunction v. Without loss of generality we can 
assume that A G (t(5+). Clearly then v vanishes in M_. Recalling (13. 9p and 
(j3.13p one obtains 

-w'l + a^wi = Q, Cg(-1,1), w[[-l) = d, w[{l) = v' 

The problem admits a unique solution, since a is not an eigenvalue of (j3.16p . 
In view of ()3.7p , (j3.14p function vi is a solution of the problems 

(5.1) -v'l + Uvi = \vi, xEM_, vi{-{)) = wi{-l), 

(5.2) -v'i + Uvi = \vi + \iv, XGIR+, vi{+Q) = wi{l) - v' 

on each half-axis. The first of problems admits a unique solution in L2(M_) 
since A cr(5_). The solution of the second one does not generally exist. 
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According to Predholm's alternative problem ()5.2p has a solution in L2(M_|_) 
if and only if Ai = v' {+0){v' (+0) —wi{l)) . To obtain this, we multiply equa- 
tion (|5.2p by eigenfunction v and then integrate twice by parts. Solution 
vi is ambiguously determined, namely it is defined up to term cv. Subor- 
dinating it to the condition vvi dx = we fix it uniquely. Recalling 
()3.10p . p.l5|) and tt; = 0, we obtain the problem 

-u;2 + a^'t«2 = 0, ^e(-1,1), 

(5.3) 

w'^i-l) = v[{-0), w'^il) = v[{+0) + v"{+0), 



which gives us the corrector W2- 
We introduce the notation 

(5.4) Ae = A + eAi, Y^ix) -- 



v{x) + evi{x), \x\ > e, 

ewi{e^^x) + e'^W2{£~^x), |x| < e 



for the constructed approximations of eigenvalues and eigenfunctions. Such 
approximations can also be found in the case, when either A G a{S-)\a{S+) 
or A G a{S-)na{S+). 

5.2. Case of Open (^'-Barrier. We now assume that a belongs to reso- 
nant set Sij,, and w = aWa, where Wa is an eigenfunction of (13.16^ and a is 
an arbitrary constant. In this case A is an eigenvalue of (I3.23P with L2(M)- 
normalized eigenfunction v. Note that all eigenvalues of this problem are 
simple. The first condition (I3.1ip yields a = ■ Then the second one 

also holds since v{+0) = 0ii,(a)f(— 0) = Wo,{l) = aWai^)- view 

of Remark 13.21 constant a is zero if and only if A G cr{S^) D a{S+). 

Problem (j3.2ip admits solution wi, because p.22p holds. This solution 
can be represented as wi = wl + aiWa, where wl is a partial solution of 
the problem and ai is an arbitrary constant. This constant will be found 
below, but first we shall construct vi. 

Function vi satisfies equation (j3.7p outside of the origin and 

Vi{+0)-e^{a)vi{-0)=gi 

with gi = wi{l) — 9y;/{a)'Wi{—l) — v'{+0) — 9xi/{a)v' (—0), which follows from 
equalities (j3.14p . Although ai is as yet unknown, constant gi is uniquely 
determined. Indeed, 

wi{l) - e^{a)wi{-l) = u'i(l) - 9^{a)wl{-l)+ 

+ ai{Wa{l) - e^{a)Wa{-l)) = wl{l) - 6^ia)wU-l). 
Next we employ (13.10p and (13.15P to derive the problem 

. . (-w'^ + a^{Ow2 = {X-U{0))w, eG(-l,l), 

\ w',{-l)=v[{-0)-v"{-0), w',{l) = v[{+0) + v"{+0). 
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Its solvability condition can be written as 

(5.6) 9^ {a)v[ (+0) - v[{-Q) = hi 

with hi = {w{-l))-^{X-U{Q)) j\w'^di-9^,{a)v"{+Q)-v"{-Q). There- 
fore vi must be a solution of the problem 




-v'{ + Uvi = \vi + \iv, a;GM\{0}, 
vi{+d) - 9^[a)vi{-Q) = gi, e^{a)v[{+0) - v[{-0) = hi. 



Free parameter Ai in the right-hand side of equation (j5.7p enables us to solve 
the problem. In light of Fredholm's alternative, (j5.7p admits a solution if 
and only if Ai = giv'{—0) — hiv{—0). For definiteness sake, the solution is 
subject to additional condition vvi dx = 0. 

Given ui, we can now compute constant ai. From the first condition in 
(I3l^ we deduce oi = (VF'„(-l))-i (vi(-O) - v'{-0) - u't(-l)). It is easy 
to check directly that the second condition in (|3.14p also holds. 

Finally in the case of open J'-barrier one obtains such approximations 
for the eigenvalue and the eigenfunction of the perturbed problem 

Ae = \ + e\i, 

(5.8) {v{x) + evi{x), \x\>£, 

~ \vSFT)^"(^"'^)+^^i(^"'^) + ^'^2(e-^x), |x| <£. 

Here W2 is an arbitrary solution of (j5.5p . Recall that condition (j5.6p ensures 
its existence. The choice of 02 in representation W2 = W2 + a2Wa is of no 
importance since we do not look for corrector V2- 

6. Justification of Asymptotic Expansions 

As shown in Theorem 12.11 for some singular potential a^s (x) there is 
finite number N{a,'^) > of eigenvalues A|(a,^'), converging to —00 as 
e — > 0. Other eigenvalues remain bounded as e — > 0. We shall show that 
these eigenvalues converge to the eigenvalues of TC{a,^). 

6.1. Convergence Theorem. Let {Agjggj be a sequence of eigenvalues 
of 7^e(a,^) and assume that {yi;}e&i is a sequence of the corresponding 
-L2(M)-normalized eigenfunctions. Here X is an infinite subset of (0, 1) for 
which is an accumulation point. 

Theorem 6.1. // Ag ^ A and ^ v in L2(M) weakly as I 3 e ^ 0, 
then A is an eigenvalue of TC{a,'^) with the corresponding eigenfunction v. 
Furthermore, converges to v in the L2(M.)-norm. 

We have divided the proof into a sequence of lemmas. 

Lemma 6.1. Under the assumptions of Theorem \ 6.1l for every 7 > 
sequence y^ converges to v in VF|(M \ (—7,7)) weakly as well as in the 
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norm of C^{M. \ (—7,7)). Moreover, v is a solution of the Schrddinger 
equation 

(6.1) - v" + Uv = \v 
on each of half- axes. 

Proof. Let be the set of test functions if £ C^{M) such that (p{x) = 
for X G (—7,7)- We conclude from (13. ip that 

(6.2) / y'^ipdx= / {U-Xe)ye<fdx 
Jm. Jr 

for all if G Ai-y and e < 7 since the support of short-range potential ^'^ 
lies in (—7,7). The right-hand side of (j6.2p has a limit as e ^ by the 
assumptions, thus the left-hand side also converges for all if £ A4^. We 
deduce then that ^ u in W2^(M \ (—7, 7)) weakly and 

/ v"ipdx = {U — X) vipdx for all ip G 
Jr Jm 

It follows herefrom that u is a solution of (jG.ip on M \ (—7,7), and so on 
whole half- lines in view of the arbitrariness of 7. It remains to note that 
the imbedding theorem implies convergence of y^ in C^(M \ (—7,7)). □ 

Lemma 6.2. Suppose as in Theorem \6.1\ that ^ A and y^ ^ v in L2(M) 
weakly as 1 3 e ^ 0. Then ye(e) ^^(+0), ^^(e) f'(-|-0), ye(— e) 
v{-0) and y'^{-e) v'{-Q) aslBe^O. 

Proof Let C be the C~(M \ {0})-function such that C{x) = for x < 
and X > 2, and also ("(x) = 1 for x G (0,1). We hereafter denote the 
characteristic function of set K by xk ■ Let us introduce sequence Ce (x) = 
X{£,oo){x)C{x). Multiplying both equations (|3.ip . (j6.ip by Q and integrating 
by parts yield 

00 00 
y'M = -J yU'dx + /(A. - U)y,Cdx, 



00 00 
v'{e) = - j v'C dx + j\x- U)vC dx. 



Here we take into account that Ce(e + 0) = 1 and C'eix) = for x G (e, 1). 
According to Lemma fe. li the right-hand sides of the equalities have the same 
limit as e ^ 0, then ^/^(e) v'{+0). In a similar way, using cut-function 
Ce{—x), we can prove that y'^{—e) v'{—0). 

Now let us consider the cut-function rj G C°°(R\{0}) such that ri{x) = 
for X < and x > 2, and also r/(x) = x for x G (0, 1). We set rj^{x) = 
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X(e,oo)(^)^(2;)- 111 view of equalities r/e(e + 0) = e and C" (^) = for x G 
(e, 1), and employing p.ip . ()6.ip we derive 



= ev'(e) — / vq" dx + {U — X)vr]dx 



1 e 

Because ^/^(e) t''(+0), we conclude that the right-hand sides of these 
equalities converge to the same limit, and so ye(e) v{+Q). The same 
proof works for the sequence e). □ 

We are now in a position to describe the behavior of eigenfunctions 
in a neighborhood of the origin. Let w and z be solutions of the Cauchy 
problems on [—1, 1] 

(6.3) -w" + a^{i)w = {), u;(-l) = 1, = 0; 

(6.4) -z" + a^(e)z = 0, z(-l)=0, z'{-l) =v'{-^). 

Lemma 6.3. If ^ X and ^ v in L2(IR) weakly as T 9 e — > 0, then 

(6.5) \\e-he{ei) - e-^ye{-e)w{i) - 2(e)||ci([_i,i]) ^ 0. 

Proof. Set tUe(C) = e~^ye(e^) — £~^y£{—£)w{S,) — z{0- Upon substituting 
^ = e~^x we see that equation (13. ip may be rewritten as 

Employing (|6.3p and (|6.4p we deduce that Wg, is a solution of the Cauchy 
problem 

lu;,(-l)=0, ^x;^(-l) = y^(-e) - z;'(-0) 



(6.6) 



with fe{i) = e{Xe - U{ei))ye{ei). For every E L2(-l,l) there exists 
a unique solution We G VF2(~1'1)- Moreover this solution satisfies the 
estimate 

(6.7) Ik.llw-K-i,!) < C + \y'A-e) - v'{-^)\) 

with constant C being independent of e. Next the inequality 

/I 1 /"^ 1 1 

vli^i) di = - I y1{x) dx < - ||ye||L2(iR) = - 
-1 £ J -e ^ ^ 

implies e^'^/'^\\yeieC)\\L2{-i,i) < c. Hence ||/£||l2(-i,i) < ce^^'^, because 
is bounded. Therefore in light of Lemma 16.21 the right-hand side of ([67 
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tends to zero. The proof is complete by using the imbedding theorem 
PF|(-l,l)cCi([-l,l]). □ 

Lemma 6.4. // Ag — > A and ye ^ v in L2(M) weakly as I 3 s ^ 0, then 
ys ^ V in L2(M). 

Proof. First we show that y^ is bounded on M uniformly with respect to e. 
Applying (j6.5p and Lemma [6. 2 1 we see at once that sequence y^ is uniformly 
bounded on [— Indeed 

(6.8) max < ce + \z{s^^x)\e + \ye{—e)\ \w{e^^x)\ < ci. 

x£[~e,e] 

Let = R \ (— e,e). Multiplying the equation (j3.ip by Xns{x)ye{x) and 
integrating by parts yield 

J y'edx = Xe J yldx- j Uyldx + y'^{-e)ye{-e)-y'^{e)ye{e), 

All terms in the right-hand side are uniformly bounded with respect to e. 
In particular, boundedness of the term including potential U follows from 
the WKB-asymptotics of y^ for large x [29| p. 55]. Hence sequence y^ is 
bounded in W2{^e)-: and so in C{Q.e)- This and (|6.8p yield max||/£(x)| < c 

with constant c being independent of e. 

Fix 7 > 0. We see from Lemma |6. II that for small e the L2(l^'y)-norm of 
difference — u is smaller than 7. Thus 

he - vWhiCm < he - v\\l2{^-,) + he " v\\l^^_^^^) 

< (1 + 2 max \ye{x) — v(x)|)7 < C7 

with constant C being independent of e. The proof is complete by recalling 
that 7 can be made as small as is wished. □ 

Proof of Theorem \6.1\ We conclude from Lemmas 16. H 16.41 that t; is a solu- 
tion of the equation 

-v" + Uv = Xv, xeR\ {0} 

and ||f ||l2(ir) — 1- We are left with the task of showing that v satisfies 
appropriate coupling conditions at the origin. In light of lemma 16.31 we 
thus deduce 

(6.9) e'^ {yeis) - ye{-e)w{l)) ^ z{l), y',{e) - e-^y,{-e)w' {I) ^ z'{l) 

as e — > 0. According to Lemma [6.2( we have ye(e) — ye{—s)w{l) 0, 
ys{—s)w'{l) 0. Therefore 

(6.10) v{+0) - v{-0)w{l) = 0, 

(6.11) v{-0)w'{l) =0. 
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First let assume that v{—0) ^ 0. Then w'{l) = 0. By ()6.3p w is a non- 
trivial solution of the boundary value problem —w"+a^w = 0, w'{—l) = 0, 
w'{l) = 0. Hence w is an eigenfunction of operator 7^ and coupling con- 
stant a belongs to resonant set Sij/. Furthermore 6xi/{a) = w{l) since 
w{—l) = 1 by construction. Thus (16.10p coincides with coupling condition 
(f3T20]l . We obtain from the second formula in (lOj) that v'{+0) = z'{l), 
and so z is a solution of the boundary value problem 

(6.12) - z" + a^z = 0, z'{-l) = v'{-0), z' (l) = v' (+0) , 

which coincides with (I3.2ip . Then solvability condition ()3.22p must hold, 
which is the second coupling condition for v. Thus v is an eigenfunction of 
T-L{a, ^) with the eigenvalue A. Here A is a limit of sequence A^. 

Now assume v{—0) = 0. Employing (|6.10p yields u(+0) = 0. If q T,^, 
then the proof is complete. In the case, when a belongs to Sij,, function 
V satisfies condition (j3.20p and w is an eigenfunction of 7^. Then there 
exists solution z of (|6.12p . consequently condition (j3.22p holds. This case 
is described in Remark 13. 2[ □ 

Corollary 6.1. Suppose that eigenvalue A|,(a,^') of TCs{(x,^) is bounded 
from below. Then A|(a, ^) has a finite limit as I B e ^ and this limit is 
a point of the spectrum ofTC{a, 

Proof. Suppose, contrary to our claim, that 

/X* = lim A|(a, < lim A|,(a, ^) = fi*. 

The constants fi^, fi* are finite since A|(a,^') is a bounded function. Re- 
call that A|(a,^') is a continuous function of e € (0,1). Then for each 
A G [fi^,iJ,*] there exists a subsequence of eigenvalues Xg, e £l, converging 
to A. For instance, set 2 can be chosen as a sequence of the roots of equation 
A|(a, ^') = A with respect to e. Sequence {ys}e(^i of the corresponding nor- 
malized eigenfunctions contains a weakly convergent subsequence. By the 
theorem A is an eigenvalue of Ti.{a,'^). Therefore interval [fi^,fi*] belongs 
to spectrum a{7i{a, ^')), but this is impossible unless = /x*. □ 

Corollary 6.2. For each eigenvalue A of Tl{a,^) with multiplicity s there 
exist exactly s eigenvalues A|(a, ^) of Ti.£{a,^) converging to X as e ^ 0. 

Proof. Operator TC{a, ^) has simple or double eigenvalues. Let us assume 
that s = 1, but A|.(a, ^) — > A and A|_,_]^(a, ^') — > A for some k. Then there 
exist two sequences {y|(x; a, ^')}egj, {y|_(_]^(x; a, ^')}£gj of the eigenfunc- 
tions, which converge in L2(M) to vectors of the form e^'^v. This contradicts 
the fact that y|(x;a,^') and yf._^_^{x;a,^) are orthogonal in -L2(M) for all 
£ € 2. Similar arguments apply to the case s = 2 (see [21j for more de- 
tails). □ 
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6.2. Quasimodes of 'He(a, ^'). We now show that each point of cr('H(a, ^')) 
is a hmit of the eigenvalues of Ti.e^a,'^). Let i? be a self-adjoint operator 
in Hilbert space H with a domain T){B). A pair (/i, u) G M x T>{B) with 
\\u\\h = 1 is called a quasimode of the operator B with an accuracy up to 
p > if \\Bu — hu\\h < P- 

Lemma 6.5. Suppose that the spectrum of B is discrete and simple. If 
(fj-ju) is a quasimode of B with accuracy to p > 0, then interval [fi — p, p + p] 
contains an eigenvalue A of B. Furthermore, if segment [p — t, p + t] 
contains only this eigenvalue of B, then \\u — v\\h < 2t~^p, where v is an 
eigenfunction of B for the eigenvalue A and \\v\\h = 1- |44l I57j 

Inasmuch as we didn't construct the asymptotics for the case of double 
eigenvalue A in the previous section, we assume that cf{S-) H = 

which ensures that the spectrum of TC{a, ^) is simple. 

Let us construct the quasimodes of ^£(0,^). Now fix some eigenvalue 
A of 7^(0, with eigenfunction v, HwllLaCK) = 1- For each A and v we have 
obtained the formal asymptotic approximations A^, defined by either 
()5.4p or (jS.Sp depending upon a and The task is to justify the asymp- 
totics. Notice that in this subsection we don't distinguish the resonant and 
non-resonant cases. By construction A^ and satisfy equalities 

-y" + {U{x) - Ae)Ye = e^Ri{£,x) for \x\ > e, 
(6.13) -y/ + {U{x) + ae^2^'(e~^x) - Ke)Ye = eR2ie,x) for |x| < e, 
[Ye].=±s = e^rfis), [y;],=±e = e^rfie) 

with functions Rj, rf^ being uniformly bounded with respect to its argu- 
ments. Function Y^ doesn't belong to the domain of '^^(q,^'), because it 
has jump discontinuities at points x = ite. But we can construct function 
with the following properties 

• Ce is a smooth function out of points x = ite and differs from zero 
only for e < |x| < 1; 

• [Ce]x=±e = -rfie) and [Cs]x=±e = -rf (e); 

• max ([^^(x)! + |Ce(a;)| + |Ce'(3^)l) ^ c with constant c being inde- 

e<|x|<l 

pendent of e, 

which eliminates this defect. In fact, function 1^ + e^C^ is continuously 
differentiable at x = ite and belongs to ^^{Tieia,^)). We set = yy^ + 
e^Ce|l2^2\R)(-^ + e^Ce) arid substitute into (|6.13p instead of YJ.. Then 
the orders of smallness of right-hand sides in (j6.13p don't change since 
lly^ + e^C£||L2(K) ~^ 1 e — > 0. Therefore pair (A^, T^) is a quasimode of 
TCeia,"^) with accuracy to e. 
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Lemma 6.6. For each A G (T('H(a, ^')) there is an eigenvalue A^(a,^') of 
7is{a,'^) such that Aj(a;, — > A. Moreover, 

(6.14) |A5(a,^') - A| < cie, \\y'j{- ;a,^) - v\\^^^^^ < C2e, 

where y| and v are the corresponding normalized eigenf unctions. 

Proof. Let (A^, T^) be the quasimode of 7is{a, ^) corresponding to the limit 
eigenvalue A and eigenfunction v. According to Lemma 16.51 there exists 
number j such that Xj{a, ^) — < ci£, from which the first inequality 

in ()6.14p follows. In view of Corollaries 16.11 [6^2] index j is independent of e. 
If T is less than the distance from A to the rest of the spectrum of Tl{a, ^'), 
then for small enough e interval [A — r, A + r] contains eigenvalue Xj{a, ^) 
only. Applying again Lemma [631 yields • ; a, ^') — TgH^^^^j^^ < 2t~^ci£, 
from which the second inequality in (j6.14p immediately follows. □ 

We are now in a position to state the main result. Let {Xk{a, ^')}^]^ be 
the eigenvalues of H{a, ^') numbered by increasing and let {vk{x; a, ^')}^^ 
be the orthonormal system of eigenfunctions. Recall that = N{a, ^) is 
the number of eigenvalues of Ttei^a, ^) converging to — oo as e — > 0. 

Theorem 6.2. Assume (a, ^') G M x P and a{S-) n cr(S'+) = 0. Then for 
every natural k we have 

(6.15) |A|+7v(a> ^) - ^fe(a> ^)| < cie, 
(6-16) ||y|+;v(-;a,^) " • ; «> ^) IL2(k) < 
with constants ci, C2 being independent of e. 

Proof. Theorem 16. H Corollary 16.21 and Lemma 16.61 may be summarized 
by saying that eigenvalue A|_,_^(a, \I') tends to Afc(a,^') as e ^ 0. Thus 
estimates (j6.15p . (|6.16p are more precise versions of (|6.14p . □ 
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